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Abstract 

Using the method of Tracy and Widom we rederive the correlation functions for 
/3 = 1 Hermitian and real asymmetric ensembles of N x N matrices with A'' odd. 

1 Introduction 

The starting point for many results concerning the spectral theory of random matrices is the 
derivation of a determinental or Pfaffian form for the correlation functions of the eigenvalues. 
Perhaps the most well-studied ensembles are Hermitian ensembles having joint probability 
density function (JPDF) of the form 



where Z a normalizing constant, w : M — > [0, oo) is a weight function and A(A) is the N x N 
Vandermonde determinant in the coordinates of A. The parameter /3 is called the inverse 
temperature parameter due to its physical interpretation in the study statistics of log-gases, 
another area where joint densities of the form (1.1) arise. When w is Gaussian and /3 — 1,2 
or 4, H. is the JPDF of the ensemble ol N x N real symmetric (/3 = 1), complex Hermitian 
{(3 — 2) and self-dual (/3 = 4) matrices, the independent entries of which are chosen with 
Gaussian density. The details of this derivation are given in [15] and the interpretation of 
Q in the study of log-gases is given in [8] . We will denote the ensemble with joint density O 
by wf3E. 

The nth correlation function of w/JE is defined by 



where A V ^ = (Ai, . . . , A„, ^i, . . . , £,N~n) and /i is Lebesgue measure on R. That is, after 
renormalization, i?„ gives the nth marginal probability density of f2. 

When (3 = 2, Fubini's Theorem together with elementary row and column operations on 
the Vandermonde determinant in the integrand lead to the determinental formula 




n=l 



(1.1) 





i?„(A) =det [ifjv(A, ,Afc) 



n 



\j,k=i ' 
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where is the kernel of a certain operator on {w dfi) . 

FoUowing pioneering work of Dyson [6] , Mehta derived a Pfaffian^ form for the correlation 
functions of the Gaussian (3=1 and 4 Hermitian ensembles [14]. This was repeated for 
general weights by Mehta and Mahoux [13], except for the case (3 — 1 and N odd. This 
last remaining case was given by Adler, Forrester and Nagao [1]. In each of these cases, the 
Pfaffian formulation of i?„ is given by 

where here, is a 2 x 2 matrix which is the kernel of an operator on L'^{w rf/i) x LF'{w d^i). 
(The matrix Kj^ depends on (3 and has a slightly different structure depending on whether 
N is even or odd). 

Much of the interest in the spectral theory of random matrices revolves around the eigen- 
value statistics as iV ^ cxd. It is for this reason that the determinental/Pfaffian formulations 
for the correlation functions are so important: N appears as a parameter in the kernel and 
in many cases, K^j be analyzed as TV ^ oo. In contrast, the number of integrations neces- 
sary to arrive at the correlation functions via their definition, increases with — a situation 
which is not as easily handled. 

These determinental/Pfaffian forms of the correlation functions have since been derived 
by various different methods. Of particular note is the method of Tracy and Widom who, 
for all cases except P = 1, N odd, derive the determinant/Pfaffian forms of the correlation 
functions using the fact that, if A is a T x A'' and B is an A?' x T matrix, then 

det(I- AB) = det(I-BA), (1.2) 

where on the left hand side of this equation I is the T xT identity matrix, and on the right 
hand side it is the NxN identity matrix. (Equation (1.2) is more generally true; a fact which 
Tracy and Widom exploit to great advantage in a number of other situations). When (3 = 2 
this identity leads immediately to the determinental correlation functions whereas, when 
/? = 1 or /? = 4, the path from (1.2) to the Pfaffian correlation functions is more convoluted. 
The reason for this extra difficulty is that (1.2) is an identity about determinants, whereas 
what is really needed is an identity about Pfaffians. Tracy and Widom get around this 
difficulty by using the important observation that, if A is an antisymmetric square matrix 
of even size, then det A = (Pf A)^. 

The extra circumlocutions necessary in the derivation of the correlation functions can 
be eliminated by using, in place of (1.2), the fact that 

Pf(C-"^ - A^BA) _ Pf(B-T - ACA^) 

PfC-T ~ PfB-T ' ^ 

where B and C are arbitrary antisymmetric 2T x 2T and 2A^ x 2A^ matrices with nonzero 
Pfaffians, and A is an arbitrary 2N x 2T matrix. The derivation of the correlation functions 
using (1.3) in the case where N is even was carried out by Borodin and Sinclair in [4, 
Appendix A]. Equation (1.3) is Rains' Pfaffian Cauchy-Binet formula [16]; a proof is given 
in [4, Appendix B]. 

The main purpose of this note is to show how (1.3) can be used to derive the Pfaffian 
form of the correlation functions in the case when (3 = 1 and N is odd, thus completing 



-"^In fact Dyson, Mehta and many who foUowcd, expressed the correlation functions in terms of Quaternion 
determinants, but in the more recent literature these are usually expressed in terms of Pfaffians. 
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the work started by Tracy and Widom. We will also show how (1.3) can be used to derive 
the correlation functions of real asymmetric ensembles in the case when A'' is odd. For TV 
even a derivation of the correlation functions of real asjonmetric ensembles using (1.3) was 
given by Borodin and Sinclair [3]. The work here will complement the existing methods 
for deriving the correlation functions of real asymmetric ensemble given by Sommers and 
Wieczorek [20] and Mays and Forrester [9]. 

1.1 Real Asymmetric Matrices 

Ginibre's real ensemble oi N x N matrices is given by R^^^ together with a probability 
measure specified by treating the entries of the matrices as independent standard normal 
random variables [11]. This ensemble is complicated by the fact that there are two species 
of eigenvalues: real and complex conjugate pairs. Among the implications of this is that 
there is no one JPDF for the ensemble. Instead, we have several partial JPDFs indexed by 
pairs non-negative integers L and M satisfying L + 2M = N; L represents the number of 
real eigenvalues in a 'sector' and M the number of pairs of complex conjugate eigenvalues. 
These partial JPDFs are given by 



where Z is a normalizing constant (that depends on N but not L and M), A(a V /3Y (3) is 
the N X N Vandermonde determinant in the variables 



When L = N and M = 0, (1.4) was established by Ginibre in his founding treatise on real 
asymmetric ensembles [11]. The general case was established three decades later indepen- 
dently by Lehmann and Sommers [12] and Edclman[7]. 

We may generalize Ginibre's real ensemble to other ensembles of real asymmetric matri- 
ces by allowing w to be another function, and we call such an ensemble the real asymmetric 
ensemble with weight function w; here we will only consider weight functions satisfying 
w{X) = w(X) for all A e C. 

Given non-negative integers i and m with i + 2m < N, we define the £, m-correlation 
function of the real asymmetric ensemble with weight function why R^^m : x C™ [0, oo) , 




(1.4) 



ai, . . . , aL, Pi, Pi, . . . , Pm,/3m, 



and w : C ^ [0, oo) is a weight function given by 




where 




where /Ui is Lebesgue measure on R and /U2 is Lebesgue measure on C. 
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One goal of the current manuscript is to show that, when N is odd, Rg^. 
written as 

_KN{Zk,Xj') KN{Zk, Zk') 



i(x, z) can be 



^^,m(x,z) = Pf 



k,k' = l,- 



where Kn is a particular 2x2 matrix kernel acting on L'^ini + H2) x L?{p-i + ^-2), the exact 
nature of which will be explained in the sequel. A similar statement is true when N is even 
[3, 4, 19, 10], and there are other existing (and arguably more complicated) methods for the 
N odd case [20, 9]. 



2 de Bruijn's Identities 



2.1 For Hermitian Ensemble 

Given a measure z/ on M, we define 



Z'^^0= [ |A(A)|V(A). 



When dv = wdfi for a function w : M [O. oo), Z'^ ^ is the normalizing constant for the 
corresponding ensemble with weight function w. We call -^jv^^ the partition function of the 
corresponding ensemble. 

A first step in the derivation of Pfaffian or determinantal form of the correlation functions 
for such ensembles is to write Z'^^^ as a determinant or a Pfaffian. This can only be done 
when /3 = 1, 2 or 4. These identities, the first due to Andreief [2] and the second two 
to de Bruijn [5], can be formulated in our setting as follows: Suppose po(A),pi(A), . . . are 
arbitrary monic polynomials such that degp„ = n, then 



Z^2 = ^! det 



-I N-l 



Pj{\)pk{X)dv{X) 



Z]^,4 = (27V)!Pf 
and when A'' is even 



/ pjiX)p'kiX)-Pk{X)p'j{X)du{X) 



-I 2N-1 



Z^NA = Pf 



/ / PjWPkiv) sgn(77 - A) dv{ri)dv{X) 



-I N-l 



j,k=0 



We define the N x N matrix ^ by 



PjWPkiv) sgn(?? - A) dv{'q)dv{X) 



N-l 



j,k=0 



(2.1) 



Thus, when N is even, Z^ ^ = Pf UJ^ One reason de Bruijn did not provide an identity 

for Z^ ^ for odd A'^ is that the Pfaffian is only defined for antisymmetric square matrices 
with an even number of rows and columns. We may produce a de Bruijn identity in the odd 
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N case by suitably altering the matrix to appear on the right hand side of the expression 
for Z'^ ^. Specifically, for A'^ odd, we have ^ = Pf ^, where 



/ 

Jr 



- [ po{X)di^{X) 



/ pN-i{X)du{X) 
Js. 



PN-i{X)di^{X) 




■ (2.2) 



2.2 for Real Asymmetric Ensembles 

The corresponding identity for real asymmetric ensembles is given in [17]. This can be 
written down as follows: Given a measure ui on R and a measure 1^2 on C \ M which is 
invariant under complex conjugation, define v = vi + U2 and set 

^N= Yl WJIm I I |A(«V/3Y;9)|dj.f(a)rfi.f(/3). 

I,+2M=JV 

If dvi = wd^ii and dv2{l3) ~ \w\'^d^2 for some fimction w : C ^ [0, 00), then Z^ is the 
normalizing constant for the real asymmetric ensemble with weight function w. When N is 
even, 

Z'^ = Pf UJ^, 

where UJ^ is the N x N antisymmetric matrix given by 

JV-l 



N 



and when N is odd, Zf^ = Pi WJ^ where 



W 



N 



Pj {(^)Vk (7) sgn(7 - a) dvi {^)dv-i (a) -'^i j Pj {P)Pk (/?) sgn Im(/3) dP2 (/?) 

c 

/ po{X)du,{X) 
Jr 

/ pN-iiX)di'i{X) 

Jr 

- I pN-i{X)di^i{X) 



j,k=0 



Po{X) dviiX) 



3 Correlation Functions in terms of Partition Functions 

From here forward we will limit our attention to /? = 1 Hermitian ensembles and real 
asymmetric ensembles. We will drop the subscripts on all relevant quantities for these 
ensembles so that, for instance, Z'^ represents both Z'ff and Zf^; which partition function 
is being represented will be clear from context. This will allow us to treat these cases 
simultaneously at the most technical part of the proof. 

We will ultimately be interested in the case when A'' is odd, but the results in this section 
are equally valid for iV even. 
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3.1 For Hermitian Ensembles 

We first consider the case of wlF,. Let T > N he ana even integer, yi,y2, ■ ■ ■ jVt € R and 
suppose ci, C2, . . . , ct are indeterminants. We define the measure 

T 

dvW = ^Ctw{yt)d6{X - yt), 
t=i 

where S is the probabihty measure with unit mass at x = 0. We also define the measure v 
by dv = wd^. Z^^'^ /Z" is the generating function for the correlation functions of wlE; this 
is the content of the following lemma. 



Lemma 3.1. Given an integer n> we define n = {1,2, ... ,n}. 
Z'n-Vv 

= 1 + 2. 



7+1/ T f " 1 

F^ = l + Z] H j n^tO) f^"(2/i(l),•••,yt(n))• 



The proof of this lemma follows easily from the definitions of Z'^^'^ and -R„; details can 
be found in [3]. 

3.2 For Real Asymmetric Ensemble 

For the real asymmetric ensembles with weight function w we suppose U and V are even 
integers greater than and set 

u V 

dvW = X] ay,w(Xu)dd{C - Xu) + by(w{Zy)dd{C - Z„) + w{Zy)d{C - Zy)), 

u=l v=l 

where xi, . . . , xjj G K, zi, . . . ,zv G C \ K and ai, . . . ,au ,hi, . . . ,hv are indeterminants. If 
1/ = i/i + 1/2 is the measure with dvx = wdjii and dv^ = |w^pc^A*2, then Z^^" /Z" generates 
the correlation functions of the corresponding real asymmetric ensemble. 



Lemma 3.2. 



X! X! X] ] H '^uW ^t)(fe) l'-R^,m(a;u(l)) ••• )a;u(^), 2^0(1)) ••• )^0(m))) 

v.:t^U X>:m/-V 1 = 1 fe=l ' 



Z^ 

e+2m.<N -' — — — 



where by convention we will take 







u:OyNj=l t>:OyNk=l 

We remark that Lemma 3.1 follows from Lemma 3.2 by setting bi = ■ ■ ■ = by = 0- The 
proof of Lemma 3.2 is found in [3] and follows directly from the definitions of and 



4 Using de Bruijn's Identities 

From here forward we will assume that N is odd. 
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4.1 For Hermitian Ensembles 

Given a measure /t on M we define the operator e*^ : L'^{k) — > L'^{k,) by 

^yW = l I /(e)sgn(A-Od«(A). 
Using this we define a skew-symmetric bilinear form on L'^{k) given by 

Jr 

Notice that U'', as given in (2.1), can also be written as 

where, as before po,pi, . . . is a sequence of monic polynomials with degp„ = n. 

Since Z'^^'^ jZ^ generates the correlation functions, we turn now to the investigation of 
the entries of W"'"^. It will be convenient to set -pn = wpn;n = 1,2,..., and to write 
{f\g) = {f\g)^ and e = e'^ (recall /x is Lebesgue measure on M). 

{Pj\Pk)ri+u = / Pj{X)Pk{0 sgn(^ - X)d{i/ + r]){X)d{y + r?)(0 

T 

= {Pj\Pk}^ + ^^ct{pj{yt)epk{yt) - Pk{yt)epjiyt)) 
t=i 

T T 

- XI ^tCmPj{yt)pk{ym) sgn{ym - yt)- (4.1) 



t=l TO=1 

The other entries in W"'"'' are of the form 



T 

/ pj{\)d{'n + v){\)= [ pj{X)di^{X) + J2ctPj{yt) 

Jw Jr 
4.2 For Real Asymmetric Ensemble 

First wc extend the definition of the sgn function to C by specifying that sgn(z) = if 
z ^R. Next, given a measure k = ki + k;2 on C we define the operator : L'^{k) ^^(k) 

by 



i / /(Osgn(A-OdKi(A) AeM; 
— isgn5(A)sgnIm(A) AeC\] 



As in the Hermitian case, we define a skew-symmetric bilinear form on L'^{k) by 

/ {fe^g-ge^f)dK. 



As before we set e and {f\g) for e** and (/Iff)'*, where in the context of real asymmetric 
ensembles, n = fJti + 112- 



8 



Correlation Functions for (3=1 Ensembles of Matrices of Odd Size 



A superficial change of notation will bring the entries of W''+'' into a form identical to 
the Hermitian case. We setT = U + V, rename the indeterminants ai,. . . , au, 61, . . . , 6y to 
Ci, C2, . . . , ct and rename xi, . . . , xu, zi,. . . ,zv to yi, 2/2, • • ■ 7 Vt- If we also define 



dS{X - 7) 7 G R; 

d6{X - 7) + d6{X - 7) 7 e C \ ] 



d6{X - 7) = 
It follows that we may write 

T 

dvW = X! ctw{xt)d5{X - yt), 

t=i 

and, as in the Hermitian ensemble case, 

T 



{pj\Pk)r,+v = {Pj\Pk)^L + 2^ct{pj{yt)epkiyt) - Pk{yt)<^Pj{yt)) 

t=i 

T T 

- X] X] (^tCmPj{yt)pkiym) sgn{ym - yt), (4.2) 



t=l TO=1 

and 



5 The Method of Tracy and Widom 



Since we have written Z"^^^ jZ" in a uniform manner for both Hermitian ensembles and real 
asymmetric ensembles we may, at least for the moment, analyze both cases simultaneously. 
We define A to be the N 2T matrix given by 



y2crpo(yi) y2c7epo(yi) 
\f^\p\{y\) \/2crepi(yi) 

/2cIpjv- 1(2/1) ^/2c[epN -\{y\ ) 
Xr(2/i)Vci/2 

where Vr is the characteristic function of K 



A„,2t-1 
A„,2t = 



'^CtPniyt), 

^epn{yt), 



n ■■ 
n ■ 



V^Poiyr) 
V^c^Piiyr) 

/2c^Pjv-i(2/t) 




. That is, 

:0,1,...,7V-1 
:0,1,...,7V-1 



V^epoiyr) 
V^epiiyr) 

^2c^€pN -i{yT ) 

X&{yT)\/cT/2 ^ 



1,2,. ..,T; 
1,2,. ..,T, 



and 



A-N,2t-i =0, t = 1,2, . . . ,T; 

Ajv,2t = XR(yt)\/cf/2, i = l,2, ...,T. 
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and set J to be the 2T x 2T matrix, 



1 
-1 



1 
-1 



It follows that, 



AJA^ 



^^ct{pn{yt)epm{yt) - epn{yt)Pm{yt)) 



N-l 



n,m=0 



t=l 



('tx^{yt)Po{yt) ■■■ - X] (^txs.{yt)pN-i{yt) 



'YctXR{yt)PN-i{yt) 





L t=l t=l 
Next, we define the 2T x 2T matrix 



E 



yJCtCu 



sgn(2/„ -yt) 











SO that, 



E2t-i,2u-i = sgn(y„ - 2/t); u = 1, 2, . . . , T, 



and all other entries of E are equal to 0. An easy computation reveals 

T T 

(AEA"^)„,„ = ^ ^ c„ctp„(y„)p„(yt) sgn(y„ -t/t); n, ro = 0, 1, . . . , TV - 1, 

u=\ t=l 

and all entries in the last row and column of (AEA^) are 0. That is. 



AEA' 



2 ^ ^ CuCtPm{yu)Pn{yt)^yt,yu) 



«=i t=i 



-, JV-l 



- n.m— 



It follows that 



y^r,+-^ =W'' + AJA^ - AEA"^. 
And, if we define the AT + 1 x TV + 1 matrices C""^ = W and B = -J + E, then 
Z-^+n ^ pf w''+'' _ Pf (C-"^ - ABA"^) 



Z-^ PfW^ 



Pf c- 



(5. 
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We are now in position to use the identity 

Pf (C-T - ABA"^) _ Pf (B-T - A^CA) 



Pf C^T 

and an easy calculation shows that 



PfB-T 



(5.2) 



B-T = -J - E' 



where E' is the 2T x 2T matrix, 

E' = 





„ \/CtC~ 



sgn(2/„ - yt) 



t,u=l 



It is also easily verified that Pf B ^ = (—1)^. Thus, 

= (-l)'^Pf (- J-E'- A'^CA) =Pf (J + E' + A'^CA). 
We now compute the entries of A^CA: 

Af-l AT-l 

(A"^CA)2„-l,2t-l = '2'\fc^ ^ ^ Pn{yu)Cn,mPm{yt); U,t = 1,2, ... ,T; 



n—0 m— 
N-1 N-1 



N-1 



(A"^CA)2„-l,2t = ^Vc^ ^ ^ Pn{yu)Cn,mePm{yt) + y/c^XRiVt) X] Pn{yu)Cn,N; 

n=0 m=0 n=0 

u,t= 1,2,..., T; 



N-l N-1 



N-1 



(aTca)2„ ,2t-l = 2a/c„C( ^ ^ ePn{yu)Cn,mPm{yt) + a/c^Xr(2/m) X/ ^ 

N,mPm (yt); 

n=0 m=0 m=0 

u,t= 1,2,..., T; 

N-1 N-1 
(A"^CA)2„,2t = y/CuCtXRiVt) ^ C„,jvePn(yu) + y/c^XKiVu) X] ^N,mePm.{yt); 



n=0 
N-1 N-1 



m=0 



+ 2Vc„Ct ^ ^ ep„(y„)C„,TOepm(yt); u, i = 1, 2, . . . , T. 

n=0 m=0 

If we define 

JV-l iV-l 

DSN{y,y') = 2 X! X! Pniy)Cn,mPm{y'); 



n—0 m— 
Af-1 Af-1 



N-1 



SN{y, 2/') = 2 X] X] Pn{y)Cn,mePm{y') + X^iv') ^ Pn{y)Cn,N; 
n—0 m— n—0 

SNi{y,y') = 2 X! X] ^Pn{y)Cn,mepm{y') 

N-1 

+ X] c„,jv(xR(y')ePn(y) - XR(y)epn(y')). 



n—0 m— 



n=0 
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then 



A^CA = 



DSN{yu,yt) SN{yu,yt) 
-SN{yt,yu) SNl{yu,yt) 



(5.3) 



u,t=l 



and, if we define the 2x2 matrix kernel Kn : C x C — > C^^^ by 

KN{y,y') 



DSN{y,y') SN{y,y') 
-SN{y, y') SNl{y, y') + \ sgn{y' - y) 



and the 2T x 2T matrix K by 

K= [^/o;IctKN{yu,yt)\[^^^-^ 
then from (5.1), (5.2), (5.3) we see that 



Pf (J + K) . 



(5.4) 



6 Recovering the Correlation Functions 



The final step in the derivation of the correlation functions is to expand Pf(J + K) using 
the identity 

T 

Pf(J+K) = l + ^ ^ PfKt, (6.1) 

where Kt is the 2n x 2n minor of K given by 

= [VCtO) Ci(fe) Kn {yt(j) , yt(k) )] ■ 
A Proof of this identity can be found in [21] or [18]. 

6.1 For Hermitian Ensembles 

Applying (6.1) to (5.4) and using the fact that 

Pf Kt,t = I H ct(„) I Pf [ifjv(2/tO),2/t(fc))]2fc=i ' 

m=l 



we find that 



~Z^ 



1 + 



Yl Y \ n '^tM f [-^Jv(2/tO),2/t(fc))]2fc=i- 



This together with Lemma 3.1 implies that 

Rn{y) = P{[KN{yi,yk)]lk=^. 
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6.2 For Real Asymmetric Ensembles 

Given u: iy U and \3 : m y V, we define u V : i + m / U + V by 



(uVt.)(i) = / "(-^'^ 



3<^ 



li £ + m = n, then each t : n y T_ can be written as uV for some u : £ y U_ and : m F . 
(This does not preclude the possibility that t = u or t = o). Using (6.1) we find 



Pf(J + K)= E E E PfKuv», 



e<u,m.<v ~ ~ — — 



where we take the Pfafhan of an empty matrix {i.e. when £ and m equal 0) to be 1. It is 
straightforward, if technical, to show that 



pfK„vo=] n «"0) n f 

i=i fc=i 



Kn (a^uO) ) a;u(i' ) ) (a;uO) , 2o (fe') ) 



j,j'=i,...,e 

fc,fe' = l,...,m 



Thus, from (5.4), 



E E E n^u^n^" 



o{k) 



0<« + m<T ~ 



k=l 



Pi 



,-f''jv(2;o(fe),a;u(fe')) -fl'jv(^B(fe), 2^B(fe')). 



This result and Lemma 3.2 imply that 

^<!,m(x,z) = Pf 



(6.2) 



KN{Zk,Xj') KN{Zk,Zk') 



j,3'=i,...,e 

fc,fc' = l, . . . ,m 



7 Simplifying the Matrix Kernel 

On first inspection, the kernel is dependent on {po,P2, ■ ■ ■ ,Pn-i)- In fact, can be 
shown to be independent of this family of moni polynomials — though a wise choice may 
simplify the representation of K^. In this section we suppose the existence of a complete 
monic family of polynomials q — {qo, gi, . . . , qn-i) such that 

{Q2n\q2m-l)i' = - {q2m-l\Q2n) u = 5n,m.rn n = 0, 1, . . . , J, 

where J = (A'' — l)/2 — 1. Here r = (n, r2, . . . , rj) are referred to the normalizations of q. 
We also define 



Sn= qn{a)di^{a), 
Jm 

Clearly q, r and s are dependent on w. 



n = 0,l,...,N -1. 
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Using these definitions, W becomes 



' 















So 


















Sl 









ri 









S2 






-n 












S3 













rj 





SjV-3 










-rj 








SN-2 











•• 











SN-1 


-So 




-S2 


-S3 • • 


• -SN-3 


-SN-2 


-SjV-1 






An easy computation reveals, 



Z- = PiW = SN-illrj. 

In order to present the entries of Kn in their most simplified form, we need to invert 
W. Recall that C = (W)""^. Exphcitly, C is given by 





1 



Sl 



1 




-So 





1 

ri 



S3 



1 

n 




-32 



roSjv-1 r-QSjv-i nsjv-i nsjv-i 




rj 

SN-2 



1 

rj 


-SAr-3 



TjSn-1 VjSn-I 





-Sl 

rosjv-i 
So 

rosjv-i 

-S3 

nsiv-i 

S2 

riSN-1 



— SiV-2 

rjSN-1 

SjV-3 

rjSN-1 


-1 



SN-1 



SN-1 





At the entry level, 

C2j,2m+1 = — C2j+l,2m = 
C2j,N-l = — Cjv-l,2j = 

and 



■^2j,2m 



■^2j+l,2m+l 



S2j+1 
rjSN-l ' 



, c 



2j+l,JV-l 



SjV-1 



^JV-l,2j+l 



m = 0, 1. 



i,m = 0, 1,..., J; 

S2j 



TiSN-l 



0,1,..., J; 
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It follows that 



DS^iyy) = 2j2'^Ml2j+iiy')-'l2j+i{y)q2j{y') 



3=0 



I ^QN-ijy) ^ S2j+iq2j{y') - S2jq2j+i{y') ^ 2qN-i{y') ^ S2jq2j+i{y) - S2j+iq2j{y) 
^ Vj SN-1 ^ rj 

J 



^ 2qN-i{y) S2j+ieq2j{y') - S2jeq2j+i{y') 

3=0 ■> 

2eqN-l{y') S2j(I2j + l(JJ} - ■'^2j + L(l2j{y) ^ gjV-i(jy)A.'.(?yO 

Sn-1 rj Sn-1 

3=0 ■' 



and 



J 



SNl{y,y') = 2j2- 



^(i2j{jj')(<i2j+i(j/) - ((}2j+i{y)i'(i2j{y') 



3=0 



2 ( i\ - i ( /y ) \^ g2j+i»]2j(,'yO - ■'^2j<q2j+iiy') 



Sn-1 . „ - J 

3=0 ■' 

^ 2eqN-i{y') S2jeq2j+i{y) - S2j+ieq2j{y) ^ eqN -i{y)xAy') " ^qN-i{y')x^{y) 
Sn-1 ^ Tj Sn-1 
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